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Abstract 


The objective of this study is to investigate the location of the gas-liquid interface at various temperatures in a polymer electrolyte membrane 
fuel cell under non-isothermal conditions. A mathematical model, coupled with the electrochemical process, two-phase flows, species transfer, and 
heat transfer is employed. A finite volume-based CFD approach is applied to investigate the species transport behavior in a fuel cell. The effects of 
two model parameters, namely cell temperature (T,.\,) and humidification temperature (T), on the gas-liquid interface and cell performance are 
presented. Simulation results indicate that variations of these two parameters influence the location of the gas—liquid interface, the cell performance, 
and the distribution of liquid water saturation. At lower cell temperatures, the gas—liquid interface moves toward the inlet port of the channel when 
the humidification temperature is greater than the cell temperature. Therefore, the cell performance decreases as the liquid water clogs the passage 
for the transport of oxygen. Furthermore, these two factors are closely related to the membrane temperature distribution. Obvious variations in 


magnitude are seen at a cell temperature of 323 K and a humidification temperature of 343 K. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Proton exchange membrane (PEM) fuel cells are currently 
regarded as promising energy-conversion devices for the 21st 
century because of their high efficiency, low emission, and quick 
start-up capability, making them suitable as power supplies of 
mobile, stationary, and portable devices. However, cell perfor- 
mance is determined by several factors, including the geometry 
and morphology of transport components, as well as the operat- 
ing conditions such as pressure, temperature, and humidification 
of the reactant gases [1-3]. 

The transport behavior in a fuel cell involves fluid flow, 
thermal transfer, and mass transport processes, along with the 
electrochemical reaction. Furthermore, water and heat manage- 
ment are critical to the proper operation of PEM fuel cells. A 
complex two-phase transport process prevails due to the pro- 
duction of water by oxygen reduction in the cathode catalyst 
layer. Excess liquid water may block the gas pore paths for 
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mass transport through the porous medium, substantially reduc- 
ing the reactant transfer rate to the reaction sites and thereby 
degrading cell performance. However, the formation of liquid 
water depends on the saturation vapor pressure, which is strongly 
dependent on the temperature. Therefore, the temperature factor 
is an inevitable consideration in water management investiga- 
tions because the phase change of water such as condensation 
and/or evaporation closely relates to the corresponding satura- 
tion pressure [4,5]. 

Several models have been developed to demonstrate the influ- 
ence of liquid water on water management. Wang et al. [6] first 
studied two-phase flow and liquid water transport on the cathode 
side of a PEM fuel cell based on the multiphase mixture model 
(M? model) originally developed by Wang and Cheng. The 
model encompassed both single- and two-phase regimes cor- 
responding to low and high current densities and was capable of 
predicting the transition between these two regimes. The half cell 
including the gas channel, gas diffusion layer, and catalyst layer 
on the cathode side of a PEM fuel cell was considered and an 
expression for calculating the gas—liquid interface was presented 
in the article. Later, this model was further developed to include 
a complete cell by You and Liu [7]. Pasaogullari and Wang 
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Nomenclature 

A specific electrochemically active area 
a water activity 

C molar concentration (mol m~?) 


Q 


quadratic drag factor 
constant-pressure heat capacity (Jkg~! K7!) 
capillary diffusion coefficient 
binary diffusion coefficient (m2 s7!) 
Faraday constant (C mol`!) 
gravitational acceleration (m s72) 
enthalpy (Jkg~!) 

current density (A m~?) 

transfer current density (A m7?) 
permeability (m7) 

thermal conductivity (W m`! K7!) 
molecular weight (kg mol~!) 
capillary pressure 

pressure (atm) 

universal gas constant (J mol`! K7!) 
rate constant 

source term 

liquid water saturation 

temperature (K) 

cell temperature (K) 

humidification temperature (K) 
fluid velocity (m s7!) 

mass fraction 

molar fraction 


sempana TRSTE mD YG 


Subscripts 

a anode 

c cathode 

con condensation 

e electron 

eva evaporation 

eff effective value 

g gas phase 

l liquid phase 

p proton 

sat saturation value 

w water 

Superscripts 

ref reference value 

Greek 

a transfer coefficient 

E porosity 

(o) phase potential (V) 

n overpotential (V) 

K proton conductivity (S m7!) 
Àk mobility of phase k 

u dynamic viscosity (N s m~?) 
v kinematic viscosity (m? s7!) 
Ac contact angle 


p density (kg cm7?) 

Pk kinetic density (kg m7?) 

o electrical conductivity (S m7 ty 

Keff ionic conductivity of the membrane (Q7! m7!) 


T tortuosity of the porous medium 
E stoichiometric flow ratio 
y membrane water content 
¢ surface tension (N m7!) 


[8] developed a model to explore the two-phase flow physics 
in the cathode gas diffusion layer. The simulations revealed 
that flooding of the porous cathode reduced the rate of oxygen 
transport to the cathode catalyst layer. Furthermore, they indi- 
cated that the humidification level and the flow rate of reactant 
streams are key parameters controlling PEM fuel cell perfor- 
mance, two-phase flow, and transport characteristics. Meng and 
Wang [9] improved their three-dimensional model based on 
recent visualization experiments, and more accurately investi- 
gated the two-phase behavior under different gas utilizations. 
Although energy transport was ignored, their model was able 
to predict liquid water flooding dynamics. Luo et al. [10] pre- 
sented a three-dimensional, two-phase isothermal model based 
on the M? formulation for investigating the condensation and/or 
evaporation interface under the condition of low-humidity inlet 
reactant gases, and they successfully predicted the important 
feature of the dry-wet-dry transition (DET) in the PEM fuel 
cell. 

Zhang et al. [11] studied liquid water transport and its removal 
from the gas diffusion layer and gas channel of PEM fuel cells 
both experimentally and theoretically. The distribution of liquid 
water on the gas diffusion layer surface and inside the gas chan- 
nel was observed in situ. They characterized the formation and 
emergence of liquid water from the gas diffusion layer surface 
and identified two modes of liquid water removal from the gas 
diffusion layer surface. Ju et al. [12] presented the first work in 
which a mathematical model of the energy equation was sys- 
tematically derived from thermodynamics. The thermal model 
was further coupled with the electrochemical and mass trans- 
port to deal with the heat and water management in PEM fuel 
cells. Later, Wang and Wang [13] expanded the M? model to per- 
mit investigation of the interaction between the two-phase flow 
and thermal transport due to non-isothermal effects. The results 
revealed that vapor-phase diffusion enhances water removal and 
provides a new mechanism for heat removal through a phase 
change process. This new heat removal mechanism is similar to 
the heat pipe effect. Shimpalee and Dutta [14] presented a math- 
ematical model to predict the temperature distribution inside a 
straight-channel PEM fuel cell and to study the effect of heat 
produced by the electrochemical reactions on fuel cell perfor- 
mance. Source terms for transport equations, heat generation, 
and a phase change model were presented to facilitate their 
modeling work. The results showed that the temperature profiles 
depend on the generation of heat by electrochemical reactions 
and phase changes of water inside the fuel cell. Further, the 
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cell performance is not merely influenced by the inlet humidity 
conditions, cell voltage, and membrane thickness, but also by 
the temperature rise within the fuel cell. Natarajan and Nguyen 
[15] developed a pseudo-three-dimensional model by extend- 
ing a two-dimensional isothermal model for conventional flow 
fields. Their results indicated that water removal by evaporation 
increases at higher temperatures and at higher stoichiometric 
flow rates, which results in better cell performance. However, 
if the water produced by the electrochemical reaction cannot 
be effectively removed from the cathode side, then flooding 
of the electrode will occur, reducing access of the reactants 
to the catalyst layer. Yang et al. [16] reported on the mechan- 
ics of liquid water transport; starting from emergence on the 
gas diffusion layer surface, droplets grow and move toward the 
mainstream flow in the gas channels. Furthermore, the oversat- 
urated water vapor causes more liquid film formation on the 
hydrophilic channel walls and the consequent channel clogging 
affects cell performance. Liu et al. [17] studied experimentally 
the two-phase flow of reactants and products as well as the water 
flooding situation in the gas flow channel of a cathode in a PEM 
fuel cell with three flow fields. Their study addressed the effects 
of flow field, cell temperature, cathode gas flow rate, and oper- 
ation time on water build-up and cell performance. The results 
indicated that liquid water forms in columns that accumulate 
in the cathode flow channels and clogs the porous medium, 
limiting mass transfer, thereby resulting in degraded cell per- 
formance. Further, the amount of water in the flow channels at 
high temperature is much less than that at low temperature. 
The above literature review clearly indicates that water and 
thermal management are critical issues in investigations of 
PEM fuel cells. However, discussions concerning the effect of 
the gas-liquid interface on the cell performance and transport 
process are rarely seen. Therefore, in this study, we report the for- 
mation and influence of the gas—liquid interface along the flow 
channel direction at various cell temperatures and humidification 
temperatures using a non-isothermal, three-dimensional, multi- 
component, two-phase model of a PEM fuel cell. Additionally, 
the temperature distribution in the cell domain membrane and 
the distributions of temperature and liquid water saturation in 
the cross-section of the cathode gas diffusion layer in the inlet 
region have been investigated. The developed model has been 
solved by means of commercial computational fluid dynamics 
software based on the finite volume approach [18-22]. 


2. Mathematical model 


In this study, the multiphase mixture formulation (M? model) 
has been adopted as it is particularly suitable for two-phase flow 
modeling in PEM fuel cells [6,8—10]. It regards the multiphase 
system as a chemical mixture and the multiphase flow is then 
described in terms of a mass-averaged mixture velocity and dif- 
fusive flux. The physical model employed in this study consists 
of nine parts, namely bipolar plates, gas flow channels, gas dif- 
fusion layers, and catalyst layers on both the anode and cathode 
sides, and a polymer membrane that is sandwiched between 
them, as illustrated in Fig. 1. Symmetry is assumed, and a single 
straight flow channel is therefore considered herein. 
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Anode 
Channel 
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Anode CL ys 
y 
Cathode CL 


Cathode GDL % 


Y Cathode 
Channel 


Cathode BP 


Fig. 1. Schematic structure of PEMFC and computational domain. 


Utilizing the M? formulation for two-phase transport, the 
proposed model incorporates the following assumptions. 


1. The system operates in a steady state. 

2. The gas mixtures at the anode and cathode are ideal 
gases. 

3. The mixture flow field is laminar as its velocity is low under 
normal operating conditions. 

4. The properties of the porous medium are isotropic and homo- 
geneous. 

5. The mass source/sink is negligible in the continuity equation 
[23]. 

6. The thermal conductivities of the cell components are con- 
stant. 


2.1. Governing equations 


The non-isothermal, three-dimensional, multi-component, 
and two-phase model of a PEM fuel cell includes five nonlinear 
coupled conservation equations of mass, momentum, energy, 
species, and charge, which are described as follows. 


2.1.1. Mass conservation equation 

Since the gaseous and liquid water are present simultaneously 
in the control volume, the conservation equation of mass for a 
multiphase mixture is: 


V - (€eprpu) = 0 (1) 


where p represents the density of the mixture and is defined as 
the volume-weighted average of the phase mass concentration 


C.-I. Lee, H.-S. Chu / Journal of Power Sources 171 (2007) 718-727 721 


in the two-phase flow [6]. When liquid water is present, the 
effective porosity is given by £eff = e(1 — s). 


2.1.2. Momentum conservation equation 

The general form of the Navier-Stokes equation is used with 
source terms that describe the drag forces in the porous medium. 
The equation is: 


V - (peep) = —EeffV P + V - (€ctt err Vit) + EeffOkg + Sm(2) 


where [Weg is the effective viscosity of the mixture, px the kinetic 
density, and Sm is the sum of the Darcy and the Forchheimer drag 
forces, 


Sm = —(Fharcy or Frorch) (3) 


which are given by: 


2 > 
es E 
Darcy drag force = FDarey = a a 
F ESCEP > > 
Forchheimer drag force = Frorch = — FAR B 


JK 


The parameters Cp and K represent the quadratic drag factor and 
the permeability. 


2.1.3. Species conservation equation 
The species conservation equation for the gas mixture is 


V > (eeft Ck) = V - (Dreff Y Ck) + Se (6) 


where k represents the chemical species, including hydrogen, 
oxygen, nitrogen, and water. Dk eff = De” represents the effec- 
tive diffusion coefficient of the kth component of the fuel reactant 
[24]. The exponent t on the porosity € is the tortuosity of the 
porous medium. The source term S. defines the production or 
consumption of the kth species in the gas phase and is given 
by: 


H? : -——— j 7 
2 2FC,"* ( ) 
02: -= (8) 
n-ar" 
H20: Lj (9) 
Dee” 
The terms ja and je are defined as: 
C 1 
. 4 ref | “Ho (aa F/RT) 
ja = Ajo (3) fe a n ETET. (10) 
2 
Co. 1 
>. __ 4 -ref 2 (@aF/RT)n _ 
je = Ajo (=) fee n FET] (11) 
2 


i 1.417(1 — s) — 2.12(1 — s} + 1.263(1 — s}, 
J(s) = 


1.4175 — 2.125? + 1.26353, 


The above conservation equations of mass, momentum, and 
species are derived on the basis of the M? model. The con- 
stitutive relationships of mixture parameter and variables are all 
dependent on the liquid saturation, defined as the ratio of the 
liquid volume to the pore volume: 


Vi 


S= (12) 

Vpore 
Density : p = pis + pg(1 — 5) (13) 
Concentration : pC = piCis + PgCg(1 — s) (14) 
Effective porosity : £eff = e(1 — s) (15) 
Kinetic density : px = pjA1(s) + pgAg(s) (16) 

e 
Effective viscosity : ef = PIs + Pal s) (17) 
(krı/vi) + (krg/Vg) 

i p ky = 8, liquid phase 

Relative permeability : 3 (18) 
krg =(1—s)°, gas phase 
k 

l P ae ac a Pr re 

Relative mobility : ky/V + krg/Vg 
Ag(s) = 1 — à(s), gas phase 
(19) 


Furthermore, during the operation of a fuel cell, liquid water 
forms if the partial pressure of water vapor exceeds the satura- 
tion vapor pressure. The liquid water thus formed may occupy 
the pores and thereby prevent the diffusion of fuel, causing 
mass transport overpotential in the porous medium. Hence, 
the effect of liquid water is taken into account. Additionally, 
capillary forces dominate the transport of liquid water on the 
hydrophilic surfaces because the pores in the porous medium 
are extremely small. Therefore, the generalized Richards 
equation, developed by Wang and Beckermann [25,26] to elu- 
cidate two-phase flow transport in capillary porous media, is 
applied: 


> Kk,k — 
V -(eanpitn(s)) = V: (spys - Eke E) 4 5 


ky Vg + krg vI 
(20) 


where D° and P* represent the capillary diffusion coefficient and 
the capillary pressure, respectively: 


Kkrikrg(dP® /ds) 


D° = 21) 
kri Vg + krgvi 
1/2 

P° = ç cos a. (“=) J6) (22) 


and where J(s) is the Leverett function, which takes the following 
form [27,28]: 


Be < 90° 
6. > 909 (23) 
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In Eq. (22), the contact angle, 0c, of the gas diffusion layer 
is dependent upon the hydrophilic (0° <8. <90°) or hydropho- 
bic (90° <6, < 180°) nature of this layer, and varies with the 
Teflon content. We assume here that the gas diffusion layer is 
a hydrophilic medium. Further, the surface tension, ¢, for the 
liquid water-air system is taken as 0.0625 N m7! [27]. 

Sı is a simplified switch function between condensation 
and/or evaporation of liquid water under these non-equilibrium 
conditions [29]. When the partial pressure of water vapor 
exceeds the saturation pressure of water, liquid water may form 
and occupy the pores in the porous medium. Conversely, the 
liquid water will evaporate if the partial pressure of water vapor 
is less than the saturation pressure of water: 

EeffXw 


Mircon 7 OwP = Psat), ifxwP > Peat 


Sı = (24) 


Teva€ettSP\(XwP — Puyo),  ifxwP < Psat 


where rcon and Feya are the condensation and evaporation rate 
constants, respectively; xw the molar fraction of water vapor, 
and Psat is the saturation pressure of water, which varies with 
temperature [30]. 


log iq Psat = —2.1794 + 0.029537 — 9.1837 
x 107>T? + 1.4454 x 107773 (25) 


2.1.4. Energy equation 

The heat generation sources in a PEM fuel cell account 
for the irreversible heat and entropic heat that is generated 
by electrochemical reactions, Joule heating that arises from 
proton/electronic resistance, and the latent heat of water conden- 
sation and/or evaporation. A generalized energy conservation 
equation is: 


n n—-l 


ket VT + S\N pDijVwjhi 
i=l j=1 


V - (EeffPCpu T) = V - 


|i -il 
-S+ 


+ ST (26) 


where keff represents the effective thermal conductivity. On the 
right-hand side of the equation, the first two terms represent the 
conduction energy and the reactant enthalpy flux; the third and 
fourth terms represent electrical-related thermal effects, and the 
last term is a source term, associated with the phase change. 


2.1.5. Charge conservation 

In a fuel cell, the potential gradient effect causes electrons 
and protons to move along individual paths. Solid phases control 
the movement of electrons. Electron transport generally occurs 
only in the bipolar plates, the diffusion layers, and the catalyst 
layers. However, ionomer phases control the motion of protons, 
which occurs in the catalyst layer and the electrolyte membrane. 
Potential fields in these two media are described as follows. 


For electrons: 


V - (GettV be) + So. = 0 (27) 


—ja atanode 
Ste = . (28) 
Jc atcathode 


For protons: 
V - (keiVGp) + Sp, = 0 (29) 


—ja atanode 
p=, (30) 
je atcathode 


where Ceff, Keff» Pe, Pp, Sp, and Sg denote electron conductivity, 
proton conductivity, electronic phase potential, electrolyte phase 
potential, and consumption rates of charge and product in the 
electrochemical reaction in the catalyst layer, respectively. 

The membrane conductivity is strongly related to the tem- 
perature and the water content wy. It is defined as the ratio of the 
number of water molecules to the number of charge sites [30]: 


1 1 
Kett(T) = (051394 — 0.00326) exp 1268 (x5 — 7) | (31) 


The water content of the membrane surface depends on the activ- 
ity of the water vapor, which also depends on the partial pressure 
of the water. Therefore, the empirical relationship between them 
can be applied: 


0.043 + 17.8la — 39.85a* + 36.00, O0<a<1 
v= (32) 
14+ 1.4(a — 1), l<a<3 
XyP 
a= (33) 
Prat 


2.2. Boundary conditions 


Boundary conditions are necessary and crucial for solving 
the above equations. They describe the operating conditions as 
well as the model geometry characteristic of the PEM fuel cell. 
The most important ones are as follows. 


Inlet boundaries: the fuel and oxidant flow rates along the 
flow channel can be described by a stoichiometric flow ratio, 
€, which is defined as the amount of reactant in the cham- 
ber gas feed divided by the amount that is required by the 
electrochemical reaction [20]. The inlet molar concentrations 
are determined by the inlet pressure, the temperature, and the 
humidity, according to the ideal gas law. 

Outlet boundaries: fully developed flow is applied. 

Walls: Neumann conditions and no-slip conditions are applied. 
Symmetric boundaries: mass flux or momentum flux has zero 
gradients. 

Electronic phase potential boundaries: fixed total cell overpo- 
tential at the outer boundary of the cathode is specified. 


ge =0  atanode bipolar plate 


Qe = Nell at cathode bipolar plate 
be _ 9 
dy 


(34) 
elsewhere 
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Table 1 

Geometrical and operating parameters 

Quantity Value Sources 

Gas channel depth/width 0.762/0.762 mm 

Shoulder width 0.381 mm 

GDL thickness 0.3mm 

Catalyst layer thickness 0.01 mm 

Membrane thickness 0.03 mm 

Fuel cell height/length 1.534/71.12 mm 

Anode/cathode pressure 1/1 atm 

Stoichiometry, at &,/&, at 1.5/3 [32] 
1.0 A cm7? 

Relative humidity of 100/100% 
anode/cathode inlet 

Porosity of diffusion and 0.4/0.4 [33] 
catalyst layers 

Porosity of membrane 0.28 [33] 

Permeability of diffusion 1.76 x 107!1/1.76 x 107}! m? [34] 
and catalyst layers 

Permeability of 1.8 x 10718 m? [34] 
membrane 

Transfer coefficient at 0.5/1.5 [34] 
anode/cathode 

Condensation rate 100s~! [24] 
constant 

Evaporation rate constant 100 atm—! s7! [24] 

Tortuosity of the diffusion 1.5 [34] 
and catalyst layers 

Tortuosity of the 3 [34] 
membrane 

Surface tension (¢) 0.0625 Nm7! [27] 


2.3. Numerical procedures 


The governing equations with their related boundary con- 
ditions were solved using a commercial code based on the 
SIMPLE algorithm for convection-diffusion problems. Details 
of the numerical solution procedure and the code are presented 
elsewhere [4,21,22]. The geometrical and operating parameters 
of the PEM fuel cell are presented in Table 1. In the simulation, 
a uniform grid distribution is used to calculate the complex elec- 
trochemical reaction and physical phenomena in the fuel cell. 
The numbers of grid nodes used for the model were 17 in the x- 
direction, 51 in the y-direction, and 71 in the z-direction. A grid 
dependence test was carried out to make sure that the solution 
obtained was independent of the grid density used. As a conver- 
gence criterion it was imposed that the normalized residual for 
each model variable is smaller than 1074 [31]. 


3. Results and discussion 


In this study, the non-equilibrium conditions of partial pres- 
sure deviation from saturation pressure have been considered. 
Our quest to delineate the effects of the gas—liquid interface 
on a PEM fuel cell has involved the following two parts. First, 
as a preliminary study, the cell temperature and humidification 
temperature were set as being equal. Second, we considered 
a scenario whereby cell temperatures were altered at a fixed 
humidification temperature. For all of the calculations carried 
out in this paper, the fuel flows were co-flows and inlet stoichio- 


F —O— Th=Tcell= 323 K 7 
F —O— Th=Tcell= 333 K 1 
—— Th=Tcell= 343 K ] 


Voltage ( V ) 


j Hili 


Current Density ( A/cm? ) 


Fig. 2. Polarization curves at various cell temperatures with equal humidification 
temperature. 


metric ratios of 1.5 and 3 were used for the anode and cathode 
sides based on a reference current density of 1 Acm~*. Further- 
more, fully humidified hydrogen and oxygen were fed to the 
anode and cathode inlets, respectively. 

Fig. 2 shows the polarization curves of the cell at equal cell 
and humidification temperatures of 323, 333, and 343 K. The 
results reveal that the cell performance improves as the tem- 
perature is increased because a higher temperature results in 
higher catalytic activity and a higher capacity for water removal 
by evaporation. Hence, increasing the temperature is helpful in 
reducing the level of flooding. Moreover, a higher cell tempera- 
ture increases the membrane conductivity and mass diffusivity, 
and decreases the mass transport resistance. Fig. 3 presents the 
effect of temperature on the location of the interface at which 
liquid water begins to condense along the cathode flow chan- 
nel at a cell voltage of 0.7 V. Indeed, the interface is defined 
as the location where the liquid water begins to condense, and 
so the saturation value behind this interface is greater than zero 
and this gives rise to the two-phase flow region. The horizontal 
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Fig. 3. Effect of cell temperature on the location of the interface where liquid 
water begins to condense along the flow channel at a cell operating voltage of 
0.7 V. 
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Fig. 4. Polarization curves for various cell temperatures at: (a) Th =323 K and 
(b) Th =343 K. 


dotted line indicates the interface between the flow channel and 
the gas diffusion layer of the cathode. The results reveal that the 
gas-liquid interface moves toward the flow channel inlet as the 
temperature is decreased, because the formation of liquid water 
depends on the saturation pressure, which is strongly dependent 
on the temperature. Liquid water is more easily and quickly 
condensed at a lower temperature. Additionally, the formation 
of liquid water may block pore paths for mass transport through 
the porous diffusion layer to the catalyst layer, thereby reducing 
cell performance. An increase in temperature increases the satu- 
ration pressure of water vapor, in turn increasing the evaporation 
capacity of the gas stream. However, although not shown here, 
a high temperature may result in the membrane drying out. 

In the second part of this study, the effect of cell temperature 
was investigated for a fixed gas humidification temperature. As 
the stoichiometric flow ratio and the humidification temperature 
do not change in this phase of simulation, the reactant flow rates 
were kept constant. Polarization curves at Th =323 and 343 K 


at various cell temperatures are plotted in Fig. 4. The curves 
in Fig. 4(a) reveal that the cell performance improves as the 
cell temperature is increased from 323 to 333 K, but deterio- 
rates as the temperature is further increased from 333 to 343 K. 
This result probably follows from the fact that the humidifica- 
tion temperature is less than the cell temperature, and the inlet 
fuel gases are therefore unsaturated. Therefore, water is almost 
entirely present in vapor form. Accordingly, the mass fraction of 
water vapor increases along the channel because of evaporation 
of liquid water from the catalyst layer. Hence, a high tempera- 
ture reduces the water content of the membrane, dehydrating it, 
and reducing its ionic conductivity. Fig. 4(b) indicates the pos- 
itive effect of cell temperature on cell performance. Since the 
humidification temperature exceeds the cell temperature, liquid 
water forms when the inlet gases enter the channel. Some of the 
liquid water keeps the membrane moist and improves its ionic 
conductivity. Accordingly, the cell performance improves as the 
cell temperature is increased. 

Some experimental investigations [35-37] have demon- 
strated that the best working conditions for the single cell 
are those under which the humidification temperature slightly 
exceeds the cell temperature. The objective of this study has been 
to investigate the impact of the formation of liquid water on cell 
performance. Therefore, the case in which the humidification 
temperature exceeds the cell temperature has been considered. 
Fig. 5(a and b) show the location of the gas-liquid interface 
along the flow channel at a cell operating voltage of 0.7 V and 
a humidification temperature of 343 K as the cell temperature is 
varied. The results reveal that the gas-liquid interface gradually 
moves to the gas inlet as the cell temperature is decreased, which 
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Fig. 5. Effects of cell temperature on the location of the interface where liquid 
water begins to condense along the flow channel at a cell operating voltage of 
0.7 V and a humidification temperature of 343 K. (a) Z=0-7.112 cm and (b) 
Z=0-0.08 cm. 
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may be ascribed to the following two reasons. (1) Reducing 
the cell temperature reduces the saturation pressure, increasing 
the amount of liquid water generated, shifting the gas-liquid 
interface closer to the gas channel inlet. (2) Some of the water 
produced in the cathode catalyst layer by the electrochemical 
reaction is evaporated because the cell temperature is higher 
and, therefore, the gas—liquid interface moves closer to the cat- 
alyst layer. Fig. 5(b) shows an enlarged portion of the curves 
along the channel from 0 to 0.08 cm. In conclusion, reducing 
the cell temperature can reposition the gas-liquid interface and 
cause liquid water to appear, which is detrimental to cell per- 
formance, since this water may occupy the pores in the porous 
medium, reducing the amount of fuel gases that can reach the 
catalyst layer. 

Fig. 6 displays the liquid water saturation field in the gas 
flow channel and the gas diffusion layer on the cathode side at 
various cell temperatures. Each of the figures reveals that liq- 
uid water saturation increases along the flow channel, because 
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Fig. 6. Liquid water saturation field in the cathode gas channel and diffusion 
layer along the flow channel at a cell voltage of 0.7 V and a humidification 
temperature of 343 K. (a) Teei = 323 K, (b) Teen = 333 K, and (c) Teen = 343 K. 


water is generated in the cathode catalyst layer by an electro- 
chemical reaction. Therefore, liquid water may first appear near 
the interface between the membrane and the cathode catalyst 
layer close to the channel outlet. Consequently, the liquid water 
saturation level in the catalyst layer is higher than that in the 
gas diffusion layer on the cathode side. Liquid water is trans- 
ported from the catalyst layer towards the gas diffusion layer 
only by capillary action; when the liquid water has reached the 
interface between the flow channel and the gas diffusion layer, 
it is transported along the channel by the drag force arising from 
the convective flow of gas. Fig. 7 compares the temperature dis- 
tributions in the membrane from 323 to 343 K. Each of these 
figures reveals that the temperature distributions in the mem- 
brane gradually decline along the flow channel, and that the 
distribution is symmetric about the z-center line. Since the inlet 
gas is fully humidified, the rate of the reaction and the temper- 
ature are higher at the inlet. The temperature variation involves 
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Fig. 7. Temperature contours in the membrane at a cell voltage of 0.7 V and a 
humidification temperature of 343 K. (a) Tce =323 K, (b) Teen = 333 K, and (c) 
Teen = 343 K. 
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irreversible heat, entropic heat, Joule heating, and latent heat. 
Additionally, the degree of temperature variation in the mem- 
brane is greater when there is a larger difference between the 
cell temperature and the humidification temperature. 

Figs. 8 and 9 show plots of the distributions of liquid water 
saturation and temperature in the X—Y section of the cathode gas 
diffusion layer in the inlet region at a cell voltage of 0.7 V and 
a humidification temperature of 343 K at various cell temper- 
atures. Fig. 8 reveals that the liquid water saturation increases 
under the land and that the amount of liquid water decreases as 
the cell temperature is increased. In the gas diffusion layer, two 
transport mechanisms operate—gas-phase diffusion as a result 
of the concentration gradient from the channel to the land and 
liquid water transport by capillary forces from the land to the 
channel. The liquid water cannot be discharged by the land, so 
the liquid water saturation on the land is higher than that in the 
channel. Additionally, the amount of liquid water decreases as 
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Fig. 8. Liquid water saturation distributions in a cross-section of the cathode gas 
diffusion layer in the inlet region at a cell voltage of 0.7 V and a humidification 
temperature of 343 K. (a) Tece =323 K, (b) Teen = 333 K, and (c) Teen =343 K. 


(a) Channel 


(b) Channel 


(c) Channel 


X 


Fig. 9. Temperature distribution in a cross-section of the cathode gas diffusion 
layer in the inlet region at a cell voltage of 0.7 V and a humidification temperature 
of 343 K. (a) Teen =323 K, (b) Teen =333 K, and (c) Teen =343 K. 


the cell temperature is increased because the saturation pressure 
increases. Fig. 9 shows that the temperature gradually decreases 
from the channel to the land. This result has two explanations. (1) 
Gas-phase diffusion occurs from the channel area to the land; (2) 
heat produced in the cathode catalyst layer by the electrochem- 
ical reaction causes the temperature gradient. Additionally, the 
degree of temperature variation is greater when there is a larger 
difference between the cell temperature and the humidification 
temperature. 


4. Conclusions 


A three-dimensional, non-isothermal, multi-component, and 
two-phase mathematical model has been developed in the frame- 
work of the computational fluid dynamics code. The effects of 
temperature on the location of the gas-liquid interface have been 
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investigated. According to the presented results and the analy- 
sis, when the anode and cathode humidification temperatures 
are equal to or higher than the cell temperature, the gas-liquid 
interface moves toward the flow channel inlet as the temperature 
is decreased. Since the formation of liquid water may block pore 
paths for mass transport through the porous diffusion layer to 
the catalyst layer, cell performance is reduced. An increase in 
temperature increases the saturation pressure of water, in turn 
increasing the evaporation capacity of the gas stream. The tem- 
perature distributions in the membrane gradually decline along 
the flow channel. The distribution is symmetric about the z- 
center line because the cell reaction rate and the temperature are 
higher at the inlet region. Numerical analysis of the results of 
this study has also shown that gas-phase fluid diffuses from the 
channel to the land and that the capillary-driven liquid water is 
transported in the opposite direction. Additionally, the degree 
of temperature variation is greater when there is a large dif- 
ference between the cell temperature and the humidification 
temperature. 
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